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Three-level atomic gradient echo memory (A-GEM) is a proposed candidate for efficient quantum 
storage and for linear optical quantum computation with time-bin multiplexing [1]. In this paper 
we investigate the spatial multimode properties of a A-GEM system. Using a high-speed triggered 
CCD, we demonstrate the storage of complex spatial modes and images. We also present an in- 
principle demonstration of spatial multiplexing by showing selective recall of spatial elements of a 
stored spin wave. Using our measurements, we consider the effect of diffusion within the atomic 
vapour and investigate its role in spatial decoherence. Our measurements allow us to quantify the 
spatial distortion due to both diffusion and inhomogeneous control field scattering and compare 
these to theoretical models. 



I. INTRODUCTION 

Information processing that harnesses the novel prop- 
erties of quantum mechanics, such as entanglement and 
superposition, can be profoundly different and, in some 
cases, much more powerful than its classical equivalent 
[2] . It is this promise that drives the development and im- 
plementation of revolutionary quantum communication 
technologies. Some of the most significant advances in 
quantum information processing have been made using 
quantum optics techniques. In particular, optical quan- 
tum key distribution (QKD) is already a proven tech- 
nique for the secure distribution of cryptographic keys via 
a shared quantum communication channel [3, 4]. How- 
ever, quantum states are fragile. They are vulnerable 
to decoherence and measurement processes that destroy 
information content. This makes the manipulation and 
storage of quantum information a significant physical 
challenge. 

To proceed further with optical quantum communica- 
tion, material systems will be required for the controlled 
storage and retrieval of quantum light fields. A 'quan- 
tum repeater' will, for example, be necessary to extend 
the range of quantum cryptosystems [5] . Proposed quan- 
tum repeater protocols operate by the generation, storage 
and transfer of entanglement among spatially separated 
quantum memories. These memories must be capable 
of coherently storing multiple quantum states of light 
for on-demand recall with fidelity exceeding the classi- 
cal limit [6]. In addition, quantum memories are also a 
key component of proposed linear optical quantum com- 
puters [7]. 

The demand for an optical quantum memory has 
brought forth a host of competing protocols [7]. Sig- 
nificant progress has been demonstrated in a number of 
operational systems that couple light fields with atomic 
ensembles. However, there is as yet no candidate which 
meets all the benchmarks required for real world appli- 
cations. 
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Quantum storage has been demonstrated using electro- 
magnetically induced transparency (EIT) [8, 9], atomic 
frequency combs (AFC) [10-12], and Raman schemes 
[13]. Electromagnetically induced transparency (EIT) 
has achieved efficiencies over 40% [14], been used to 
store light pulses in a solid state system for multiple sec- 
onds [15], and can preserve entanglement [16] and optical 
squeezing [17-20]. Entanglement has also been stored in 
AFC memories [21, 22], which have the convenience of 
large signal band widths (several GHz) [23] . 

In this work our quantum memory is a warm- gas gra- 
dient echo memory (GEM). GEM, which is a variant 
of the controlled reversible inhomogeneous broadening 
(CRIB) protocol [24], is a photon echo memory technique 
that uses applied electric or magnetic fields to reverse 
the time evolution of an atomic coherence in an inho- 
mogeneously broadened sample. Coherent light storage 
has been achieved in room temperature rubidium vapour 
with efficiencies as high as 87% using this method [6]. 
State-independent verification using the conditional vari- 
ance and signal transfer coefficient has shown unambigu- 
ously that the memory performs beyond the quantum 
no-cloning limit [25]. 

In this paper we examine the spatial multimode prop- 
erties of a GEM. Multimode capacity describes the num- 
ber of optical modes (spatial, temporal or frequency) 
that can be stored in a memory. A quantum memory 
with the capacity to simultaneously store orthogonal spa- 
tial modes is valuable because it allows parallel storage 
and processing of multiple signals- a single-cell device 
for multi-qubit memory. Repeaters composed of multi- 
mode memories can increase the channel bit-rate dra- 
matically by multiplexing between modes [26-28] . In ad- 
dition, quantum correlated images may themselves form 
the basis of new quantum information protocols [29, 30]. 
A quantum memory with high spatial fidelity is necessary 
for the storage of such signals. 

Previous work has shown storage of images using EIT 
[31-34] and four- wave mixing [30, 35]. We also note very 
recent work demonstrating the storage and recall of con- 
secutive images in a ^^Rb vapour using GEM [36]. In 
this work we investigate the capacity of a gradient echo 
memory to store complex spatial modes and multimode 
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FIG. 1. (a) The off resonance Raman level scheme used for A-GEM. The probe beam is detuned from the transition |1) |3) 
by A, the single photon detuning. The control beam is detuned from the transition |3) |2) by Ac. The detuning relevant 
to the Raman transition is the two photon detuning S = A — Ac. (b) The properties of the ensemble as a function of the two 
photon detuning S. (b-i) The unbroadened Raman absorption feature, the medium is most opaque to light that is resonant 
with the Raman transition ^ = 0. (b-ii) The inhomogeneously broadened Raman feature of the ensemble. The magnetic field 
makes the detuning a function of longitudinal position in the ensemble, the absorption feature is broadened because a wider 
range of frequencies are resonant with some component of the ensemble, (b-iii) The real component of the ensemble electric 
susceptibility x- The dispersion is greatest at the edges of the broadened Raman feature. 



images, and measure the deterioration of spatial fidelity 
and recall efficiency as a function of storage time. We also 
perform an in principle demonstration of in-memory spa- 
tial processing. Finally, the spatial fidelity of our quan- 
tum memory is an important diagnostic tool, by exam- 
ining the deterioration of the storage efficiency for com- 
plex spatial modes we probe the impact of diffusion in 
the vapour cell memory. 



II. THREE LEVEL GEM 

The quantum memory used in these experiments is a 
three level (A) GEM which stores optical information in 
the long-lived coherence between hyperfine ground states 
of warm ^^Rb atoms. The key advantage of three-level 
memories is that they harness the long coherence times 
between negligibly coupled ground states [37]. 

A strong classical control beam couples the probe sig- 
nal to the two ground states |1) and |2) via an excited 
state 1 3) in an off-resonance Raman configuration illus- 
trated in Fig. la. For a weak probe beam and large 
single photon detuning (A) the excited state population 
is negligible during the storage process. In this way the 
probe signal is coherently transferred to a spin coherence 
Pi2(r, t) in the ensemble. The control field needs to be on 
for read and write operations, but can remain off during 
storage. 

The Zeeman sub-level splitting of the atomic medium 
can be controlled using magnetic fields. We use two 
solenoids to produce 'read' and 'write' fields with lin- 
ear gradients r] and —r] along the signal propagation axis 
of the atomic medium. These fields inhomogeneously 
broaden the resonant frequency of the Raman transition 
producing the absorption feature shown in Fig. lb. In- 



dividual frequency components of the input are mapped 
longitudinally along the cell to produce a spatial spin 
wave corresponding to the Fourier spectrum of the input 
field envelope. 



Once the wave-packet is stored, the atomic dipoles pre- 
cess with an angular velocity proportional to their local 
resonant frequency. Over time the spin excitation ac- 
crues a spatial phase variation which makes coherent re- 
emission impossible. The signal can be retrieved by re- 
versing the system's time evolution. Switching the mag- 
netic field gradient from r] to —77 at some time r after 
storage inverts the local detuning and causes the dipoles 
to rotate in the opposite direction. At time 2r the dipoles 
realign and the atomic spin wave is once again in phase. 
A time reversed echo of the original signal will be re- 
leased in the forward direction. GEM can operate with 
recall efficiency approaching unity in the forward direc- 
tion [38, 39]. 



In the A configuration pulse emission can only take 
place if all atomic dipoles oscillate in phase while the 
Raman control beam is switched on. Leaving the con- 
trol field off suppresses the echo and the signal remains 
stored in the atomic excitation. This condition can be 
exploited to control the order in which stored pulses are 
retrieved, thus enabling resequencing of time-bin qubits 
[1]. Furthermore, the applied magnetic fields make in- 
memory spectral manipulation possible [40]. These de- 
grees of control readily extend to spatial processes and 
make A-GEM a promising candidate for spatially multi- 
plexed memories. 
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FIG. 2. Experiment schematic, described in text. PBS: polarising beam splitter, EOM: electro-optic modulator, AOM: acousto- 
optic modulator, CCD: charge coupled device (camera). 



III. EXPERIMENT 

The essential schematic of this experiment is illustrated 
in Fig. 2, it shows the configuration of the beams, mod- 
ulators, cavities and detectors that prepare our gradi- 
ent echo memory. The probe and control beams are de- 
rived from a single continuous wave Ti:Sapph laser source 
which produces 1.1 W of power at 795 nm. The laser out- 
put is tuned to the ^^Rb Dl transition (F = 2 ^ F'= 2) 
using an external rubidium reference cell. The source 
laser is then blue detuned from this transition by approx- 
imately 1.5 GHz by monitoring the reference cell fluores- 
cence. 

To produce the probe we transfer the source beam 
power into FM sidebands at ±6.8 GHz with a fibre cou- 
pled EOM. We use a ring cavity (Cavity 1 in Fig. 2) to 
select the +6.8 GHz sideband from the modulated signal. 
This sideband frequency corresponds to the ground state 
splitting of ^'^Rb. 

We operate a second ring cavity (Cavity 2 in Fig. 2) 
as a spatial mode cleaner, it transmits the probe light 
in only one spatial mode selected by cavity alignment. 
The output port of the cavity is also used to combine the 
probe and control beams. 

The memory cell is an anti-reflection (AR) coated 
pyrex cylinder 200 mm long and 25 mm diameter con- 



taining a mixture of isotopically enhanced ^^Rb and 0.5 
Torr of krypton buffer gas. The Raman transition be- 
tween hyperfine states is inhomogeneously broadened by 
magnetic coils that apply a linearly varying Zeeman shift 
along the propagation axis of the vapour cell. Current 
is switched between two solenoids of opposite pitch to 
create magnetic fields for the read and write operations 
(coils A and B in Fig. 2). A third coil (not shown) pro- 
vides an adjustable constant field offset. 

After the memory, the control beam is filtered from the 
signal by absorption in a second gas cell containing a nat- 
ural mixture of rubidium isotopes. The vapour pressure 
in the memory and filter cells is fixed by adjusting the 
temperature with electrical elements. The memory cell is 
kept at 70° C and the filter cell is kept at approximately 
140°C. 

Once the filter cell has suppressed the control beam at 
the memory output, flip mirrors can be used to switch 
the signal between three detection mechanisms. In the 
first a local oscillator beam in a heterodyne configuration 
is used to perform phase or amplitude measurements on 
the signal retrieved from the memory. In the second a 
photodiode can be used to record the temporal intensity 
of stored pulses independent of spatial mode. Lastly, 
we can send the beam to a high-speed CCD camera for 
spatial mode analysis. 
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The Grasshopper 2 CCD camera from Point Grey Re- 
search features a Sony ICX285 1.4 megapixel image sen- 
sor with resolution 1384 x 1036. The total size of the 
CCD sensor is 20 x 14 mm. Each image is the average 
of fifty 30 /is exposures triggered externally from a dig- 
ital control station. Alternating images of the echo and 
background are taken so that background signal (includ- 
ing control beam leakage) can be removed dynamically. 
The subtracted background is the mean of two no-signal 
images taken before and after the echo. The response 
of the CCD was calibrated against the response of the 
photodiode to ensure that the integrated CCD measure- 
ments gave an accurate measure of total energy in the 
pulse. 

The cavity locking, probe and control beam inten- 
sities, magnetic field configuration, heater and camera 
triggers are controlled by a script specially written in 
Lab VIEW® [41]. 

IV. RESULTS 
A. Fundamental mode efficiency 

Storage time is a critical parameter for quantum 
memory systems and it is limited by multiple decoher- 
ence mechanisms. Ground state decoherence, collisional 
broadening, scattering processes and diffusion all con- 
tribute to the relaxation of the atomic spin wave. To 
improve the achievable storage time it is vital to under- 
stand the mechanisms behind the loss of efficiency and 
fidelity. 

A spatial mode investigation gives us a new window 
to the operation of our memory. In this first section we 
compare heterodyne and CCD measurements of the re- 
cah efficiency of TEM-00 echoes. The CCD, unlike the 
heterodyne detector, has neither fine temporal resolution 
nor spatial mode sensitivity. It detects all spatial modes 
and frequencies and integrates the signal over an expo- 
sure time of 30 /is. However it allows us to collect spatial 
information and to operate with complex spatial profiles 
which are not accessible to the heterodyne system. In the 
following sections we will use the CCD to investigate the 
storage of both high-order Hermite-Gauss spatial modes 
and multimode images. 

Under certain conditions degenerate four-wave mixing 
(FWM) between probe and control fields configured as 
per Fig. la and a third conjugate field can be expected 
to cause amplification of the probe signal [42] . The phase 
matching condition of the FWM process imposes restric- 
tions on the propagation of the generated conjugate field 
so that when linearly polarised coupling and probe fields 
propagate through the medium at a non-zero crossing 
angle we observe the conjugate field as a distinct spatial 
mode at the CCD camera. For this experiment we use 
circularly polarised light and operate at a temperature 
and one-photon detuning such the conjugate field inten- 
sity and the associated probe amplification are negligible. 
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FIG. 3. Comparison of TEM-00 storage efficiency measured 
by heterodyne detection (blue) and CCD camera (green). The 
heterodyne efficiency is the area under the echo pulse com- 
pared to the area of the input. The CCD efficiency is the total 
intensity on the camera compared to the input. Inset: De- 
modulated heterodyne traces showing the shape of the photon 
echo at several recall times. 



By this method we ensure that the results presented here 
are free from efficiency or spatial mode distortion due to 
FWM. 

The efficiency with the control ffeld both on and off 
during storage has previously been characterised using 
heterodyne detection [6] , but this method will report ex- 
tra losses if components of the echo have changed mode. 
The inset in Fig. 3 shows the temporal profile of TEM- 
00 echoes as a function of storage time. For higher order 
spatial modes we used the photodiode detector (which 
is not mode sensitive) to check the temporal shape of 
the recalled pulse. The time profiles of the higher order 
mode echoes were not more distorted than the fundamen- 
tal mode traces shown in Fig. 3. 

The results of efficiency measurements made using the 
heterodyne and CCD are compared in Fig. 3. The con- 
trol beam is left off during storage to minimise control 
scattering losses (see Section IV-C). The efficiency mea- 
sured by the CCD tracks above the concurrent hetero- 
dyne efficiency and the difference increases with storage 
time. The difference is light which reaches the detector 
in a spatial mode orthogonal to the local oscillator or has 
spread outside of the local oscillator beam. The increas- 
ing gap suggests that there are processes in the memory 
that disperse spatial information during storage. 
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FIG. 4. (a) Recall efficiency of control-off (green squares) and control-on (blue circles) echoes as a function of storage time 
from the CCD. The solid red line is a prediction of the control-on echo efficiency from the additional decoherence due to control 
field scattering (Section IV-C). (b) Area of the Gaussian mode echo as a function of storage time from the same data. The 
points plotted are cP' where a is the standard deviation of the echo intensity distribution. Data is shown for both control-on 
(blue circles) and control-off (green squares) storage along with solid lines showing a linear fit to each dataset. The gradient of 
this fit corresponds to the diffusion coefficient D implied by the mode expansion (shown on plot). The red line is a theoretical 
prediction for the control-on width from simultaneous diffusion and inhomogeneous control field scattering (Section IV-C). 



DifTusion of TEM-00 



The dominant spatial effect in warm vapour memories 
is atomic transport by Brownian motion. It has been 
shown that the ballistic motion of warm atoms can co- 
herently distribute a collective excitation throughout a 
gas cell [43]. Therefore we expect to measure expansion 
of warm GEM photon echoes corresponding to diffusion 
of the atomic excitation in the memory. 



Under certain conditions the internal atomic degrees 
of freedom are essentially decoupled from the atom's ex- 
ternal motion [43], and the two may be considered sepa- 
rately. The exceptions to this principle are atom-atom 
and atom- wall collisions, field inhomogeneity, Doppler 
broadening and other velocity effects [44]. In A-GEM 
memories the field is homogeneous on the plane trans- 
verse to the solenoid axis. Displacement on this plane 
does not change the internal state of each atom, there- 
fore the atomic motion can effectively be modelled by the 
addition of a classical diffusion term pdiff = DV^ p to the 
Maxwell-Bloch equations for the collective atomic states 
where D is the diffusion coefficient of the vapour [44]. 
The expectation values of the probe field envelope 8 and 
spin coherences pi2 and pis evolve according to [38] 



d£ ign i 2 n 

Pi3 = ig^ + i^cpi2 - ^ (27 + 70 + 7c) Pl3 + ^Api3 



(1) 



P12 = ^^cPi3 - (^70 + 7c - + ^ J P12 + DV'^pi2 

Where g is the vacuum Rabi frequency of the probe 
light mode, n is the rubidium density in the cell, /cq is the 
wavenumber of the probe mode and 7, 70 and 7c are the 
excited state decay, dephasing and population exchange 
rates respectively. The first equation has been simplified 
by transforming into a frame moving at the speed of light 
along the longitudinal axis z. 

In GEM longitudinal and transverse diffusion must be 
considered separately because of the applied longitudinal 
field gradient. We can model transverse diffusion with a 
Gaussian propagator. When the atomic mean free path 
is much smaller than the radius of the bounding region, 
the solution for the ensemble atomic spin wave p(r, t) far 
from the cell boundaries is the convolution of the solu- 
tion in the absence of diffusion Pstat(r,t) with a diffusion 
propagator G(r, t) of width adiff = \/2Dt [32] 



G(r,t) = {AnDt) 
p{r,t) = jdv'G{Y 



N/2 



r',t)fttat(r',0) 



(2) 
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Where N = 2 is the number of dimensions relevant for 
transverse diffusion. The recahed field is determined by 
the locally averaged value of all displaced atoms within 
a small region. That is, by the diffused atomic coherence 
operator above. The recalled field envelope mapped from 
the diffused coherence is 



\E{r,t)f 



Pit) I dr'G{r ■ 



■r',t)E{r,0) 



(3) 



Where we have introduced a homogeneous power loss 
function P(t) which includes the profile's exponential de- 
cay due to scattering and decoherence as well as losses 
caused by longitudinal diffusion which are not in gen- 
eral exponential. For a Gaussian beam (TEM-00) with 
a uniform phase front and waist Wq this diffusion model 
causes the total observed intensity to drop like 



(4) 



Figure 4a shows the decay of the recalled energy as 
a function of storage time with and without the control 
field left on during storage. Previous work has shown 
that control-off storage is more efficient as it reduces the 
amount of spontaneous Raman scattering [6] . Our model 
accounts for a spatially inhomogeneous control field and 
predicts precisely the additional losses from the control 
field. 

Using the CCD we are now able to measure the effect of 
the control beam on the mode profile of the recalled echo. 
Defining a as the standard deviation of average beam 
profile, Fig. 4b shows the area (a^) of a 2D Gaussian fit 
to the recalled pulse image both with and without the 
control beam left on during storage. 

Convolving the TEM-00 profile with a Gaussian dif- 
fusion propagator produces a broadened Gaussian with 
width cTecho = \/<^in + <^diff)- Therefore the measured 
echo intensity profile has width cr^^ho/^ ~ VK(z)^/4 + 
(7^-ff/2 where W(z) is the beam width of the input mode. 
This allows us to infer the diffusion coefficient from the 
measured echo expansion. 



pressure at room temperature [45]. This implies that 
D ^ 31 cm^/s in our memory vapour cell. 

The expansion rate calculated with the control beam 
on (Fig. 4b) is 240 ± 18 cm^/s, a factor of eight higher 
than the expected atomic diffusion rate in the cell. It 
is clear that the control beam is causing additional dis- 
tortion of the profile during storage. We explore this 
mechanism in the following subsection, and show that 
control beam scattering provides a good explanation for 
the control-on TEM-00 expansion data. 

In contrast when the control beam is off we measure 
only half the expansion expected from diffusion. In the 
control-off data the rate of mode expansion corresponds 
to a diffusion coefficient D = 13.2 dz 1.6 cm^/s. To 
resolve this discrepancy we performed a direct, inde- 
pendent measurement of diffusion in the cell using the 
method of Gozzini and Bicchi [46, 47] which indicated a 
diffusion coefficient of D = 65 =b 10 cm^/s. This appears 
to confirm that the measured expansion of the recalled 
echo signal is smaller than the expansion expected from 
atomic diffusion. 

Echo mode expansion is only an indirect measurement 
of diffusion in the cell. Indeed, there are mechanisms 
that could cause a discrepancy between the diffusion rate 
and the measured mode expansion rate. For example, 
magnetic field inhomogeneity in the transverse dimen- 
sion could reduce rephasing efficiency at the edges of the 
diffusing spin wave. Numerical simulations show that dif- 
fusion in the presence of transverse magnetic field vari- 
ation does cause additional echo efficiency loss and re- 
duced echo mode expansion. However, measurements of 
the field variation in our setup indicate field fluctuations 
that are orders of magnitude too small to explain the ob- 
served discrepancy, so we can rule this mechanism out 
as a major factor. Another possible explanation of the 
smaller than expected beam size is a reduction of the 
beam divergence due to diffusion [48]. The beam waist, 
which lies in the cell, can expand due to diffusion leading 
to a beam with reduced divergence that when measured 
downstream will appear smaller than otherwise expected. 



The diffusion coefficient is related to the collision rate 
7coii, mean free path A and atomic velocity v by 

D = Xv/3 = ^V37coii (6) 

The collision rate is well known for a number of low 
pressure binary mixtures, therefore we can estimate theo- 
retically the diffusion coefficient within the rubidium cell. 
For low pressure rubidium in krypton buffer gas the col- 
lision rate is approximately 17 MHz/Torr of buffer gas 



C. Control field scattering 

With our imaging experiment we can see that in addi- 
tion to the overall loss of efficiency the spatial profile of 
the recalled pulses is significantly worse when the control 
has been left on during storage. Figure 5 compares the 
recalled profile of several Hermite-Gauss modes with and 
without the control beam during storage. 

A significant proportion of this distortion can be ex- 
plained by the inhomogeneity of the control beam. As 
the control field power varies across the Gaussian profile 
so does the control field two-photon scattering rate F. 
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FIG. 5. The spatial profile of several Hermite- Gauss mode in- 
puts (left) and photon echoes stored for 12 /is with the control 
beam off (middle) and on (right) during storage. The image 
intensities have been normalised, i.e. they do not show the 
overall decrease in recall efficiency. 
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for A >> 



Where 731 = 27r x 5.6 MHz is the excited state de- 
cay rate [49], A = 1.5 GHz is the one photon detuning 
and Vtc is the control field Rabi frequency. For our con- 
trol beam power (400 mW) and waist {Wc = 3 mm) the 
Rabi frequency is ftc = 72 MHz. Spatially dependent 
scattering will burn out features from the probe signal 
in regions of higher control intensity. The control beam 
has a Gaussian intensity distribution which will flatten 
the recalled pulse profile, leading to an increase in the 
observed width. 

In Fig. 4a we have plotted the TEM-00 recall efficiency 
(measured by the CCD) as a function of time for control- 
on and control-off storage. The red curve is the control 
on efficiency expected from additional exponential decay 
at the control field scattering rate T. Control field scat- 
tering accounts very precisely for the additional loss of 
efficiency, and is the dominant loss mechanism when the 
control beam is on. 

Furthermore, we modelled the impact of simultaneous 



diffusion and control field feature burning and found that 
this accounts for the broadening in the TEM-00 echo pro- 
file when the control beam is on. The red curve in Fig. 4b 
is the predicted expansion from simultaneous diffusion 
and inhomogeneous control field scattering at the rate T. 
The model featured perfectly aligned co-propagating con- 
trol and probe beams with respective widths Wp = 1.5 
mm and Wc = 3 mm. The diffusion constant used in 
this model was the rate implied by the expansion of the 
control off TEM-00 mode, D ^ 13 cm^/s, but the broad- 
ening due to diffusion in this case is negligible compared 
to the broadening from the control field. 



D. Hermite-Gauss mode efficiency 

In the absence of atomic motion and control beam 
inhomogeneity we expect all spatial profiles to behave 
identically in the memory. Simulations performed with 
three-dimensional Maxwell-Bloch equations (equation 1) 
confirm that under these conditions A-GEM is insensitive 
to spatial mode. However, in the presence of atomic mo- 
tion there are reasons to expect complex spatial modes 
to decohere more rapidly during storage than a simple 
Gaussian profile [32, 50, 51]. 

During the memory write operation the transverse spa- 
tial profile of the electric field is transferred to the atomic 
excitation. This preserves spatial phase by creating a 
spin wave that mirrors the optical field envelope. When 
atomic transport is possible, atoms may drift between 
out of phase regions of the excitation. Mixing compo- 
nents of the spin wave in this way reduces the average 
coherences pi2 and pis and prevents efficient rephasing 
of the photon echo. 

The set of Hermite-Gauss spatial modes produced by 
our second ring cavity is ready-made for the investiga- 
tion of spatial spin wave interference. We stored multi- 
ple Hermite-Gauss modes for up to 60 ps and measured 
the recall efficiency of each mode as a function of storage 
time. We found that the decay rate is correlated to the 
phase complexity of the spatial profile. Figure 6 shows 
the efficiency loss curves of four higher order Hermite- 
Gauss modes TEM-mn compared to the efficiency of the 
Gaussian mode TEM-00. The echoes were stored with 
the control beam off and the control beam alignment was 
optimised separately for each mode. The magnetic fields 
were kept consistent between experiments. Each curve is 
shown with a profile of the corresponding spatial mode 
in which darkness corresponds to field intensity and hue 
indicates the field phase as given by the phase wheel at 
top right. 

Our control station sampling rate limits trigger con- 
trol precision to 1 /is. This restricts the storage times 
for which we could take data and causes leading-edge 
clipping on some echoes. The effect of this clipping is 
to alternately reduce and increase the measured inten- 
sity. The effect is periodic, every three or four points the 
echo aligns with the optimal control window. The data 
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FIG. 6. Recall efficiency of the Her mite- Gauss TEM modes with control off during storage. Control window clipping adds 
experimental noise in short cycles. The dotted line shows a three point mean filter to remove this noise. The blue trace is an 
exponential fit to the data of the form Ae~^^'^ . The decay time r from this fitting is shown on each plot. The red curve is the 
TEM-00 decay curve shown for comparison. 



is plotted with a three point smoothed function (dashed 
blue curve) which counteracts the clipping and an expo- 
nential fit to the dataset (blue curve). The red lines are 
the TEM-00 exponential fit decay curve normalised to the 
first measured efficiency of the higher order mode echoes, 
it is included to aid comparison between the figures. 

We observe decay rates that vary with the phase com- 
plexity of the spatial mode, the more complex the mode 
the faster the efficiency loss. Additionally, the decay rate 
is correlated with the distance between out of phase peaks 



in the spatial mode of the input and output fields, con- 
sistent with the effects of spin wave diffusion under the 
model of equation 3. 

There are other possible explanations for the mode- 
decay correlation observed. Magnetic field variance off 
the beam axis may become significant as mode size in- 
creases and will further reduce recall efficiency. Longitu- 
dinal diffusion should also cause super-exponential decay 
in this storage scheme. 

The Hermite- Gauss modes are not in general self simi- 
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Storage time (/is) 

FIG. 7. The inside/outside peak height ratio of control-ofF 
TEM-20 photon echoes as a function of storage time. Diffu- 
sion between the out of phase peaks causes the intensity to 
decrease fastest at the centre of the profile. The large uncer- 
tainties are due to fluctuations in the relative heights between 
images. The red line is a diffusion model with D ^ 13 cm^/s 
and also includes inhomogeneous control field scattering dur- 
ing the read and write operations. Inset: spatial profile of 
TEM-20 before (blue) and after (green) distortion induced by 
our diffusion model. 

lar under diffusion [48] , and the changing spatial intensity 
distribution provides additional evidence that spin wave 
diffusion is the dominant source of decoherence during 
storage. The TEM-20 mode, for example, has a small 
central lobe that is out of phase with the outer peaks. We 
expect, therefore, that the central lobe will be particu- 
larly prone to diffusion of atoms from surrounding re- 
gions. If we model the decay of the central lobe peak as a 
function of time (including both diffusion and inhomoge- 
neous control field scattering during the read and write 
operations) we find that it decays faster than the outer 
peaks, as shown by the red line in Fig. 7. This is a good 
fit to the experimental data, which although noisy, also 
appears to show accelerated decay of the central lobe. 
Transverse magnetic field inhomogeneity would, in con- 
trast, lead to greater decay of the outer lobes. 

We have observed that the Hermite- Gauss TEM modes 
do not decay identically. Furthermore, diffusion is a likely 
cause of the variation in the overall decay rate of the 
modes, and can be used to explain some of the finer 
structure of the decay curves with reference to the model 
introduced in Section IV-B. 



E. Spatially selective recall 

We demonstrated spatially selective recall and storage 
by aligning the control beam onto a single side of the 
TEM- 10 spatial profile. The control beam makes possible 
the storage and re-emission of the illuminated half of the 
spatial profile. The dark-control component of the signal 



profile passes through the ensemble without interacting 
with the memory. When the atomic excitation rephases 
the illuminated side of the profile is released as a photon 
echo. 

Figure 8 shows the input and echo profiles for the case 
where we have selected the left and right halves of the 
profile. This process is a rudimentary example of a spa- 
tially selective memory operation. In this experiment 
only a single half of the profile was stored and recalled. 
Spatial multiplexing in which components of a stored pro- 
file are recalled separately using only in-memory opera- 
tions is feasible with existing beam shaping techniques. 
Spatial mode operations of this sort are the basis of pro- 
posed multimode quantum networks [52] that utilise op- 
tically entangled co-propagating spatial modes [53]. 



F. Image storage 

The final result we shall present is a demonstration of 
multimode image storage in A- GEM. We present three 
images stored for up to 8 /is and analyse the operation 
of the memory as a spatial frequency low pass filter with 
a steep frequency response. 

In the previous sections the spatial profiles stored were 
pure Hermite-Gauss TEM modes generated by misalign- 
ment of a ring cavity. In order to store an image con- 
sisting of several co-propagating modes it is necessary to 
remove the redundant mode cleaner cavity and instead 
produce the desired image by passing the probe through 
a transparency mask which shapes the transmitted pro- 
file. 

In the previous configuration the probe and control 
beams were combined on the face of the mode-cleaning 
cavity exit port. For these results we combined the 
shaped probe beam with the Gaussian control beam on 
a non-polarising beam splitter. Naturally this sacrifices 
power in both beams, but we are able to compensate with 
power from the unnecessary cavity locking beam. 

The control beam was the same Gaussian beam used in 
the previous experiments. It may be more efficient to op- 
timise the control beam shape for the input image, or to 
shape the combined beam before the memory. However, 
for practical applications of image storage in multiplexed 
quantum repeaters it is important that the memory op- 
erate without prior knowledge of the input pulse. 

The images and echoes are shown in Fig. 9. We found 
the most successful image storage with a probe beam pro- 
file half the diameter of the beam used for single mode 
storage. At this beam size the images are close to the res- 
olution limit of the CCD. The storage efficiency achieved 
is very low for both images, at 4 jas the efficiency was 
approximately 4%. The spatial frequency components 
decay rapidly, consistent with the action of the diffusion 
Gaussian as a low pass filter. The diffusion timescale is 
substantially reduced by the smaller beam. 
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2.5 mm 



1 cm 

FIG. 8. (a) Input profile of the TEM-10 mode (b) The left 
half of the profile, recalled after 6 /xs storage, (c) The right 
half of the profile, recalled after the same time. 



FIG. 9. Images stored in GEM: Bat signal (left), ^ (right). 
All images have been normalised to the same intensity scale. 
Recall efficiency ^ 4% at 4 /xs. 



V. CONCLUSIONS 

Our interest in multimode storage is driven by the 
desire for multiplexed quantum memories. These are 
essential for the construction of high bit-rate quan- 
tum repeaters for long distance quantum cryptography. 
Gradient echo memories with three-level atomic ensem- 
bles have the advantage of large spatial mode capacity 
and unique signal processing capability. Using a fast- 
triggering camera we have demonstrated the storage of 
multiple spatial modes in A- GEM, as well as the storage 
of multimode images, and investigated the factors limit- 
ing spatial fidelity. 

Brownian motion within the warm rubidium memory 
cell is the dominant cause of spatial decoherence during 
control-off storage. We describe a model of diffusion in 
the memory and show how diffusion causes the atomic 
spin wave to expand during storage. We measure this ef- 
fect as broadening of the recalled photon echo. The mode 
expansion rate we measure corresponds to a diffusion co- 
efficient of D = 13.2 cm^/s, smaller than both theoretical 
predictions an independent experimental measurement. 

Inhomogeneous control field scattering is the dominant 
cause of spatial decoherence during control-on storage. 
The two-photon scattering rate depends exponentially on 
the local control field power and therefore the control field 
profile contributes to loss of spatial fidelity over time. 
This effect is present during the read and write operations 
and needs to be considered during any spatial memory 



operations. 

We demonstrated the storage of several higher or- 
der spatial modes with complex phase profiles and ob- 
served profile dependent efficiency loss rates. Diffusion 
between regions of the atomic excitation with opposite 
phase causes interference and additional efficiency loss 
consistent with our measurements. In any warm vapour 
memory diffusion of this sort poses a fundamental limit 
on the storage time for which a given profile is topologi- 
cally stable. 

Finally, we performed an in-principle demonstration 
of in- memory spatial operations using the control field. 
Spatially dependent storage and recall is merely the most 
basic of the spatial operations that are possible using A- 
GEM. The development of fast control beam shaping and 
transverse field control will provide additional degrees of 
freedom for spatial multiplexing and processing in GEM. 
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